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We study the response of random singlet quantum critical points to local perturbations. Despite 
being insulating, these systems are dramatically affected by a local cut in the system, so that 
the overlap G = |(^s|^a)| of the groundstate wave functions with and without a cut vanishes 
algebraically in the thermodynamic limit. We analyze this Anderson orthogonality catastrophe in 
detail using a real-space renormalization group approach. We show that both the typical value 
of the overlap G and the disorder average of G a with a > 0 decay as power-laws of the system 
size. In particular, the disorder average of G a shows a “multifractal” behavior, with a non-trivial 
limit a —>• oo that is dominated by rare events. We also discuss the case of more generic local 
perturbations and generalize these results to local quantum quenches. 


I. INTRODUCTION 

One of the most basic questions about quantum dy¬ 
namics is the response of a system to a sudden change 
of the Hamiltonian in a small region in space. This local 
quantum quench can be expected to generate a superpo¬ 
sition of ground and excited states, even if the system 
was prepared in its ground state before the quench, and 
the structure of this superposition is quite complicated in 
many-particle systems. An early experimental example 
of a local quantum quench from solid-state physics is the 
excitation of a core level by an incident photon of high 
energy, which appears to the conduction electrons as a 
sudden change in the local electrostatic potential. The 
response of the Fermi sea to this change is quite interest¬ 
ing 1 2 and led to the development of important concepts 
in many-body theory such as that of an orthogonality 
catastrophe^, which is a decay with system size in the 
overlap of ground states of the initial and final Hamilto¬ 
nian. 

The orthogonality catastrophe appears in the vacuum- 
to-vacuum transition probability (also known as ground 
state fidelity or ground state overlap): what is the prob¬ 
ability that the system after the quench is in the ground 
state of the new Hamiltonian? In metallic systems, this 
leads to an overlap that decays as a power-law in system 
size L, 

G=\(* b \A> a )\(xL- ( *. ( 1 ) 

In words, the orthogonality catastrophe in a metal is 
that the ground states before and after a local poten¬ 
tial change are actually orthogonal in the thermodynamic 
limit. In experiments, the orthogonality catastrophe is 
essential for understanding the frequency dependence of 
optical absorption near the X-ray edge where such ab¬ 
sorption by a core level becomes energetically allowed, 
and is also closely related to non linear I — V character¬ 
istics in quantum dots, or the Kondo effect!^ in magnetic 
alloys. The main goal of this paper is to study the orthog¬ 
onality catastrophe that emerges in strongly disordered 
quantum systems in one dimension. 


In recent years, quantum quenches in translation- 
invariant one-dimensional systems have been studied an¬ 
alytically with considerable success as powerful methods 
from confor mal field theory**^, integrabilit}® 1 ^, and nu¬ 
meric J - 11 12 ^ are available. Systems with random impuri¬ 
ties are in general more complicated but a rare example 
of a strongly disordered interacting system that can be 
studied a nalyt ically is the random-singlet quantum crit¬ 
ical point 13 HI, whose disorder-averaged equal-time cor¬ 
relation functions scale as simple power-laws but whose 
dynamical properties are relatively complicated^. Two 
unexpected findings of the present work are that there 
is a power-law orthogonality catastrophe similar to that 
in metals, even though the system is insulating, and 
that the orthogonality catastrophe is “multifractal” in 
the sense that different powers of the disorder-averaged 
overlap G scale with nontrivially different powers of sys¬ 
tem size. 

The main quench we discuss in detail is a “cut” in the 
system that disconnects the left and right halves. This 
is a physical cut, z.e., a change in the Hamiltonian, not 
the mathematical division of the Hilbert space into two 
parts that is used to calculate entanglement. However, 
there are connections between the orthogonality catastro¬ 
phe and the disorder-averaged entanglement properties 
of the random singlet phase, which have been an active 
subjectPH The crucial point is that G keeps track of how 
many singlets are affected by the cut in the system - a 
quantity different from entanglement whose statistics can 
be accessed using real space renormalization group tech¬ 
niques. The result that the orthogonality catastrophe at 
the random-singlet critical point is similar to that in a 
metal, but that there is a difference in the multifractal 
properties, are perhaps believable in light of the fact that 
entanglement entropy is known to scale similarly to that 
of a metal (or other conformally invariant system), while 
the entanglement spectrum behaves differently 17 . 

The remainder of this Introduction reviews briefly 
some relevant recent progress on related questions. New 
tools for quantum quenches include, on the experimen¬ 
tal side, atomic systems where the absence of phonons 
and resulting long decoherence times lead to many ways 
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to generate quantum quenches and observe them while 
the dynamics remain quantum-mechanical. On the the¬ 
oretical side, progress in time- depe ndent density-matrix 
renormalization group method^ 18 * 19 * allows simulation of 
systems that are neither integrable nor conformally in¬ 
variant and hence difficult to study analytically. Strongly 
disordered interacting systems remain challenging for nu¬ 
merics but are of great current importance, partly be¬ 
cause the same real-space renormalization-group (RSRG) 
methods we apply have been important for recent 
progress in many-body localization (MBLjPnmii 

the ex¬ 
istence of localized behavior at nonzero temperature or 
energy density (see Refs. 1241251 for recent reviews). While 
the ground-state fidelity is quite a different property than 
MBL, it is hoped that some of the technical developments 
here will prove useful. 

The multifractal properties of one-particle wavefunc- 
tionsjn disordered systems have been studied for many 
yeard^EZl, including recent progress on new sets of ex¬ 
ponents appearing at edges and corner^. This problem 
is loosely connected to the random-singlet phase of spin 
chains in the following way. Essentially the same random- 
singlet phase appears both for the interacting Heisen- 
berg/XXX chain as for the XX spin chain, but only the 
latter can be mapped via the Jordan-Wigner transfor¬ 
mation onto a non-interacting one-dimensional electron 
hopping system with particle-hole symmetry in the en¬ 
ergy spectrum resulting from sublattice symmetry. Using 
this trick for the XX case, the orthogonality catastrophe 
we find can be reformulated as a statement about a differ¬ 
ent kind of multifractality that appears in the comparison 
of one-electron wavefunctions between two locally differ¬ 
ent critical Hamiltonians, and in that representation is 
closer to Anderson’s original study. We wish to point 
out that the RSRG calculations we report are expected 
to be valid even for interacting systems in the random- 
singlet phase, although the ability to check them against 
microscopic numerics is limited to the XX case, and even 
then rather nontrivial. 

The remainder of this paper is organized as follows. 
Section |TT] introduces the key features of the random- 
singlet phase and the RSRG techniques used in Sec¬ 
tion [HI] to calculate the ground-state overlap. Section [TV| 
discusses the multifractal spectrum of the orthogonality 
catastrophe using a generating-function approach. Sec¬ 
tion |V| presents numerical results on the main quantities 
of interest, and Section |Vl| discusses which aspects of our 
analysis are expected to be general to other perturbations 
and to practical quantum quenches. 


II. GROUNDSTATE OVERLAP AND STRONG 
DISORDER RENORMALIZATION GROUP 

A. Local perturbations and overlaps 

Groundstate overlaps, or fidelities, are a clear-cut way 
to characterize phases of matter. They are particularly 



FIG. 1: Left: Example of random singlet configurations on 
L = 100 sites with (|4>b), top) and without (|4>a), bottom) 
a cut in the middle of the system (dashed vertical line), for 
the same disorder configuration. Only the red singlets are 
affected by the cut. Right: The overlap G = is 

computed by counting the number of loops formed by glueing 
these configurations together. In this example, we find G = 
(|) consistent with n e — 4 singlets crossing the cut in |Ta) 
and n = 6 singlets modified by the cut (see eq §). 


useful to probe the response of a system to local per¬ 
turbations: denoting by |4/a) and |4/#) the groundstate 
wave functions without and with a local perturbation, re¬ 
spectively, how does the overlap (fidelity) G = |(T^|Ta)| 
vary as a function of system size L (considering for sim¬ 
plicity a ID system)? This question goes back to the idea 
of Anderson orthogonality catastrophe, and in the case 
where the local perturbation corresponds to a cut in the 
system, is somewhat related to the bipartite entangle¬ 
ment entropy, while at the same time being fundamen¬ 
tally different and much simpler conceptually. Study¬ 
ing the response of a system to a local perturbation is 
also very natural from the point of view of local quan¬ 
tum quenches: starting at time t = 0 from the ground- 
state |Ta) of the system without perturbation and time- 
evolving with the Hamiltonian Hb of the system with the 
local perturbation, how does the time-dependent overlap 
G{t) = A \ e lH At e -iH B t |^r A ) (known as the Loschmidt 
echo in the literature) behave at long times? In this dy¬ 
namical setting, this problem is closely related to the X- 
ray edge singularity problerrP®l, an( t the Loschmidt echo 
can be related to optical absorption spectra 29 32 l For 
some recent physical applications of the Loschmidt echo 
to cold atoms or condensed matter systems, see Refs. 33 

m 

We will mostly focus on the case of a cut in the sys¬ 
tem, and study the “static” response G = |(T^|Ta)| as a 
function of system size L - we will also discuss the case of 
more generic local perturbations and come back to local 
quantum quenches in Sec. |VIB| We will ignore all phase 
factors and implicitly focus on the magnitude of the over¬ 
lap G = |(T j b|Ta)| = (T^|Ta). To be more precise, let 
14 /a) be the groundstate of a system of size L = 2 AT, and 
14/#) be the groundstate of the same system cut in two 
halves [0, N] and [V, L = 2 N]. For gapped ground states, 
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the overlap G = remains finite in the thermo¬ 

dynamic limit. For gapless systems however, G vanishes 
as a power-law as L —>• oo, with some exponent related 
to the central charge of the corresponding conformal field 
theory^. For strongly disordered systems in ID, the av¬ 
erage of G is non-vanishing because of Anderson local¬ 
ization— - although it was pointed out recently that lo¬ 
calized systems do suffer from a “statistical” version of 
the orthogonality catastrophe 41 . In this paper, we will 
be interested in random singlet quantum critical points, 
that are examples of infinite randomness fixed points. 
Despite being insulating, they exhibit critical properties 
such as alg ebraically decaying averaged correlation func- 
tion^HEHMI, logarit hmic s caling of entanglement!®, and 
energy-length scaling 43 

ln-t-iy (2) 


length scales, and iteratively constructs the groundstate 
in terms of “random singlets”. Although we do not 
expect this procedure to be accurate initially for fi¬ 
nite disorder, we will see below that the effective dis¬ 
order strength grows under renormalization, so that the 
method is said to be asymptotically exact - i.e. is be¬ 
lieved to give exact results for universal quantities. 

Let Cl = maxf { Ji} be the largest coupling in the Hamil¬ 
tonian, and let us parametrize the couplings as /% = In y-. 

We also define a RG flow parameter T = In ^ where Uq 
is the initial value of Cl. Using the decimation rule ^ 
and ignoring a factor In 2 (irrelevant at strong disorder), 
one finds the flow equation for the distribution of the 
couplings ^ 

- = ^§^+M0)Pr*Pr(/3), ( 6 ) 


(with = 1/2 for the examples treated in this paper) 
instead of the usual quantum-critical relation E ~ L~ z . 

For concreteness, we will focus on the antiferromag¬ 
netic spin-^ random-bond Heisenberg chain 

H = Y J JiSi-S i+ 1 . ( 3 ) 

i 

More generally, our results will also apply directly to the 
anisotropic XXZ chairf® 

H = E-h ( S i S i +1 + S t S i +1 + AS i S U i) - ( 4 ) 

i 

for —1/2 < A < 1, and we will discuss generalizations to 
other random systems (namely, the transverse field Ising 


chain 4 - and 


anyomc 


chains®) in Sec. 


HID 


B. Real space renormalization group 


where P r *-Pr(/?) = / dPid/3 2 Pr(Pi)Pr(P2)S W ~ /?i ~ P 2 ) 
denotes the convolution. This equation has a remarkably 
simple solution 14 ! that is essentially an attractor to all 
initial distributions 


W) = f*- 0/T - (?) 

Most features of random singlet critical points follow di¬ 
rectly from this fixed point distribution, in particular the 
scaling © between distance and energy, which can be 
recast as VI ~ r where £ is the size of the singlets cre¬ 
ated at energy scale T. It is also possible to argue that 
even though the typical value of the spin-spin correla¬ 
tion function (§o.S r ) decays as e -c ^, its average decays 
much more slowly as 1/r 2 as it is dominated by rare 
events where the two spins belong to the same singlet. 
It is also clear from this distribution that the effective 
disorder strength is T, and is therefore increasing at low 
energy. 


Disorder is a relevant perturbation (in the renormal¬ 
ization group sense) to the pure spin-| Heisenberg chain. 
At low-energy, the disorder strength flows to infinity and 
the universal properties of the random bond Heisenberg 
chain can be accessed using a real space renormaliza¬ 
tion group (RSRG) approach®, valid at strong random¬ 
ness. The key idea of this approach is to focus on the 
strongest bond of the chain Ji. Assuming strong disor¬ 
der, this bond is typically much larger than its neighbors 
Ji Ji-i,J i+ i, so to leading order we can diagonalize 
this strong bond and form a singlet between the spins 
Si and Si+ 1 , and then deal with the rest of the chain 
perturbatively. Virtual fluctuations induce an effective 
Heisenb erg c oupling between the spins Si -1 and 5^+2 
given bjf®® 

Ji-iJi+i 

j ”- - < 5) 

with J e ff <C Jf_i, Jf, Jf+i at strong disorder. Repeat¬ 
ing this process produces singlets at increasingly long 


C. RSRG, entanglement entropy and groundstate 
overlap 

The random singlet structure of the groundstate can 
also be used to infer the scaling of more involved quanti¬ 
ties of interest like entanglement entropy or overlaps. For 
example, considering a system of length L = 2 N with N 
even and open boundary conditions, the bipartite entan¬ 
glement entropy between the right and left parts of the 
system (each of size N ) is given by S = n e In 2 where n e 
is the number of singlets crossing the entanglement cut 
in the groundstate wavefunction |T^). Using RSRG, the 
disorder average of this quantity was shown® to scale as 
n e ^ | In L (here and in the following, O denotes the 
disorder average of O), leading to the critical-like scaling 
S ~ lnL. The groundstate overlap G = (^a\^b) can 
also be evaluated using the RSRG procedure: for a given 
distribution of disorder, the groundstates |\Em) and |T#) 
without and with a cut in the middle of the system are 
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FIG. 2: (a) The singlets crossing the cut in |4>a) are obviously modified in |4>s), but other singlets can also be affected by the 
cut - in this example, the singlet between the spins Si and Si . (b) Schematic representation of the typical history dependence 
in the RG responsible for the scaling of Gt yp ~ T _ln2 / 4 . Dashed lines represent singlets of |Ta) that are modified in \^b) 
(solid black lines) because of the cut. We also show in green the average RG times ln(r*/Fj) between the creation of singlets 
at energy scales T i and Tj . 


given in terms of a collection of singlets (see left panel 

of Fig. |l]) . The overlap then reads G = (^) iVloops N ^ 2 ^ 
where Ai oops is the number of loops in the configuration 
obtained by representing singlets by semicircles and by 
gluing the mirror image of |4>a) on top of |4>#} (see left 
panel of Fig. [l]). This quantity seems obviously more 
complicated than the number of singlets n e crossing the 
cut characterizing the entanglement entropy, but we will 
show that it is nevertheless possible to understand its 
scaling analytically in many cases. 

We remark that even though we consider open bound¬ 
ary conditions and a single cut in the middle of the sys¬ 
tem, the results below can be straightforwardly general¬ 
ized to a system with periodic boundary conditions for 
4 /a > cut into two halves for 4/#: all the exponents derived 
below should then be multiplied by 2 as this “periodic” 
setup now involves two cuts instead of a single one. 


III. TYPICAL DECAY OF THE 
GROUNDSTATE OVERLAP 

We first argue that the typical value of the overlap 
defined as G typ = e lnG decays with system size as a 
power law Gt yp ^ L~ atyp with <a typ = pp. Our start¬ 
ing point will be to express the overlap in terms of sim¬ 
pler observables whose statistics can be understood us¬ 
ing the random singlet fixed point. We first note that 
if the entanglement entropy of |4/ a ) is 0 - that is if 
the number of singlets crossing the boundary n e = 0 
- then G = (4f#|4/A) = 1. In general however, G is 
not simply related to n e , although it is easy to show 

that G < Q) ne ^ 2 , where we recall that n e is even. This 
is consistent with the intuition that the wave-function 


overlap G measures how different the wave-functions are 
with and without a cut, and has therefore no reason to 
be directly related to entanglement in general. For the 
random-singlet configurations that are generated by the 
RSRG, we find that it is possible to express G as 

°-U) • (8) 

where n is the number of singlets by which |4 '!a) and |4/#) 
differ, and we recall that n e is the number of singlets 
crossing the cut in |4/a), so that the bipartite entangle¬ 
ment entropy in |4/a) reads Sa = n e ln2 , while Sr = 0 
by definition. In other words, n measures how many sin¬ 
glets are modified by the cut. Note that by definition, n 
is larger than n e since singlets crossing the cut in |4/a) 
have to be modified in 14/#), so that n — n e /2 > n e / 2. 
We emphasize that eq. (|8| does not hold for arbitrary, 
generic random-singlet wave functions |4/a) and 1 4/#): it 
is an exact result that relies on the form of the (random- 
singlet) quantum states generated by the renormalization 
process and on the specific structure of the differences be¬ 
tween 14 /a) and |4 '!b) that we will describe below. 

We also note that intuitively, G essentially behaves as 
a partition function or a correlation function - this is 
actu ally how it can be calculated in the absence of dis- 
ordei® 39 f 47 l 48 l _ an( i contrary to say, the entanglement 
entropy, it is not self-averaging, with the natural quantity 
to average being InG instead of G itself. This intuition 
is confirmed by eq. (|8| , which shows that G behaves in a 
way that is closely related to the exponential of the entan¬ 
glement entropy, which we know should be self-averaging. 
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A. Scaling of n e 

In order to understand the universal scaling of G typ , 
we need to compute the disorder average of n and n e 
for a system of size L. The average n e was computed 
in Ref. [ 44 j where it was shown that the rate of singlet 
formations across the cut in |T^) is given by 

= ( 9 ) 

where g = In ~ is the total RG time between T 0 and 
T. The average RG time between singlet formations is 
therefore given by Jt = J 0 °° dgf(g)g = 3 , so that the 
number of singlets crossing the cut for a system of length 
L reads n e ~ ~ where we have used the ran¬ 

dom singlet scaling T ~ y/L between distance and en¬ 
ergy. From this, we get the scaling of the entanglement 
entropy 6 ~ ^ In L for a system with open boundary 
conditional 


B. History dependence of the RG 


The last ingredient we need is to understand how the 
average number of singlets modified by the cut n >n e 
scales with system size. In order to do this, we need to 
understand precisely the history dependence of the renor¬ 
malization procedure, by following the RG flow in the 
configurations with (B) and without (A) the cut simul¬ 
taneously. The concept of “history dependence” of the 
RSRG was also used in Ref. EH and simply refers to the 
influence of a given decimation on future RG steps. For 
example, computing the factor 1/6 in the scaling of the 
entanglement entropy n e ~ I In L requires taking into 
account the history of the singlet formations across the 
cut precisely: after being decimated, the renormalized 
central bond is typically much weaker (by a factor e _r ) 
than the other bonds in the chain and thus has a much 
smaller probability of being decimated again. 

Let us imagine that |T^) contains n e singlets cross¬ 
ing the cut, with n e even in our setup. For the same 
disorder configuration with the cut, let us focus on, 
say, the left half of the system and denote the spins 
that were involved in those n e singlets in |T^) by 
Si, 6 2 , ■ • •, S 2 i-i, S 2 i, • • •, 6 ne , where i = l,...,n e /2 la¬ 
bels the spins by increasing distance to the cut. These 
spins will be reorganized in new singlets in |T#), perhaps 
involving other spins that were not involved in singlets 
crossing the cut |T^) (which is why n can be larger than 
n e ). Focusing on the spins 622-1 and 621, two different 
scenarios can occur in |T#}: (i) 622-1 an d 62 i can form a 
singlet, (ii) 622-1 can form a singlet with another spin Sf 
to its left, which was in turn involved in a singlet with an¬ 
other spin 6/ to its left in |T^). This spin 6/ will have to 
form a singlet with a new spin to its left in IT#}, that can 
be either 622, or another intermediate spin Sf that was 
in a singlet with a spin Sf in \^a) etc (see Fig. 2a). Im¬ 
portantly, the spins modified by the cut can occur only 


between the spins 622-1 and 622, not between 622 and 
622+1 or 622-2 and 622-1* Note in particular that 622-2 
and 622-1 or 622 and 622+1 cannot form singlets in |T#), 
whereas 622 and 622-1 can. This particular structure is 
actually why (|8| holds. 


C. Scaling of n 

We now estimate how many singlets in \^a) 
(Sf , 6/), (Sf , Sf )... between 622-1 and 622 are modified 
by the cut. We imagine running the RG until the scale 
To at which 622-1 would form a singlet crossing the cut 
in | T^). Denoting by Y 1 > Yq the scale at which 622 
would form a singlet crossing the cut in |T^), we know 
that the average RG time between those two events in 
|T^) is given by In br =3 (see above and Ref. | 44 |) . At 
the scale To in the configuration B with the cut, the spin 
622-1 is effectively at the right boundary of the left half 
of the system, and is coupled to its neighbor on the left 
by a bond of coupling strength (3 given by the probabil¬ 
ity distribution Qr 0 (P) = p^e _ ^/ r °. We now ask when 
this bond is decimated - i.e. when 622-1 forms in singlet 
in | T#). Following Ref. 44 , we construct a flow equation 
for Qr(P) with the convention that J 0 °° d/ 3 Qr(/ 3 ) = pr be 
the probability that the bond involving 622-1 to its right 
was not yet decimated at scale Y. We find that Qr(/ 3 ) 
satisfies 

+ p r(0) (Pr * Qr(P) - Qr(+), 

( 10 ) 

where the first term accounts for the change in (3 when Y 
changes, and the second term corresponds to the flow of 
the coupling when the (single) neighbor of 622-1 forms a 
singlet. This equation is readily solved and we find 

Qr(P) = Y e ~ P/F 9 (m = ln ■ with g(n) = e -/ \ 

( 11 ) 

so that pr = g(fj1). The average RG duration after which 
the bond involving 622-1 is decimated in |T#) is thus 

given by ln ^ = J 0 °° dgg(ii)ii = 1 so that typically 
Yi > T a > To where T a is the scale at which 622-1 forms 
a singlet with a spin Sf in |T#). We next consider this 
spin Sf at scale Y a in the configuration A without the 
cut. Sf is then coupled to its right neighbor by a cou¬ 
pling crossing the cut that will be decimated^ at scale 
r 1 with ln p- = 2, and to its left neighbor with a cou- 

pling distributed according to Qr a (P) = y-e -/3 / Fa . Since 
Sf is situated to the right of 622, the right bond cannot 
be decimated (otherwise it would form a singlet cross¬ 
ing the cut), and the average RG duration after which 
the bond left of Sf is decimated in |T^) is again given 

by ln ^r- = / 0 °° dp J g(g)p J = 1 . This singlet in |T^) in¬ 
volves Sf and another spin Sf, that will in turn be at 
the right boundary of the left half of the system in the 
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configuration B with the cut. By using the same argu¬ 
ment, we find that Si- will form a singlet \^b) with a spin 

Sf at a scale T c given by In = 1. However, because 

In |+ = In + In + In ^ = 1 + 1 + 1 = 3, we can 
identify T c = Ti and Sf = ^ so that the cycle stops. 
After the scale Ti at which forms a singlet crossing 
the cut in |T^), all the decimations are identical in the 
configurations A and B until the scale T 2 (given typically 

by In = 3) at which another singlet crosses the cut in 
I Va). 1 

We therefore end up with a picture of the typical his¬ 
tory dependence summarized in Fig. [2bJ where within a 
total RG time In ]+ = 6, n = 4 singlets were modified by 
the cut: n e = 2 singlets crossing the cut in as well 

as two additional singlets involving the spins Sf and Si¬ 
on the left half of the system, and similarly on the right 
half. Therefore, we find n ~ ^ | InL ~ 2n e , up to 

non-universal contributions. 


D. Scaling of Gt yp 

Gathering these different ingredients, we are ready to 
compute the average lnG = — In 2 (n — n e /2). Our re¬ 
sults then imply a power-law decay of the typical value 
of G given by 

Gtyp ~ L -ln2 / 4 . ( 12 ) 

Note that although we have focused our analysis on the 
Heisenberg <© and XXZ chains Q, the results can be 
readily generalized to other random singlet spin chains. 
For instance, let us consider the random bond transverse 
field Ising chain 


It is also instructive to compare our result to the pure 
(disorder free) case, where the overlap between the wave 
functions with and without a cut can be computed using 
Conformal Field Theory (CFTThe overlap then goes 
to zero in the thermodynamic limit as G ~ L -c / 16 , where 
c is the central charge of the corresponding CFT - with 
c = 1 for the Heisenberg chain or c = 1/2 for the Ising 
chain. Note in particular that our results for the random 
case are not obtained directly from the pure case by re¬ 
placing c by the “effective central charge” introduced in 
Ref !44] for the entanglement entropy. 

IV. MULTIFRACTAL ORTHOGONALITY 
CATASTROPHE SPECTRUM 

In the previous section, we showed that the typical 
value G typ of the wavefunction overlap G = 
suffers from an Anderson orthogonality catastrophe and 
decays as a powerlaw of the system size. It is natural 
to ask whether this orthogonality catastrophe also holds 
for the mean value G. More generally, it is interesting to 
study the behavior of G a with a > 0. In this section, we 
will show that 

~ L~^ a \ (15) 

where the multifractal exponent fi{a) depends very non- 
trivially on <a, with in particular an interesting limit 
fi{ 00 ). We emphasize that we are using the word “multi- 
fractal” in a broad sense here, meaning that the exponent 
n(a) depends non-linearly on a. 

A. a —> 0 and a —»• 00 limits 


H =* - X J i (T i (T i +1 + h i a i ’ ( 13 ) 

i 


where Ji and hi are random with In hi = In J^, corre¬ 
sponding to a random quantum critical point separating 
a paramagnetic phase and a ferromagnetic phase. This 
critical point can be conveniently interpreted as a ran¬ 
dom singlet critical point where the singlets are formed 
between Major ana fermion The only difference in our 

/ \n—n e /2 

analysis is that (|8| now reads G = ( , where 


d = y/2 can be interpreted as the quantum dimension of 
a Majorana fermion, whereas we had d = 2 for the spin- 
\ chains |3| and Q. Therefore, the typical value of the 
overlap G for the critical random Ising chain (13) scales 


G+ ng ~£- ln2 /8. (14) 

More generally, our result applies directly to the random 
SUk( 2) anyonic chains studied in Ref. 05] by replacing 
factors of 2 (for Heisenberg) or \[2 (for Ising) by the 
appropriate quantum dimension. 


Using eq. ©, we first express G a as 


G a 


X P ( n ’ n e) 

n,n e 


1 \ ot{n-n e / 2 ) 

V 


(16) 


where P(n, n e ) is the probability to have a random singlet 
configuration with n singlets modified by the cut, and n e 
singlets crossing the cut in |T^). There are two limits 
that can be understood easily. First of all, in the limit 
a 0, ~ G? vp - L-« ln2 / 4 , so that 


H(u) 


In 2 

~ a-p, 

a —>-0 4 


(17) 


using the results of the previous section. 

It is also interesting to study the opposite limit a -+ 00 . 
As a -+ 00 , the only configurations that survive in 
eq. (16) are those with n = 0 and n e = 0. That is to 
say, the limit a 00 counts the number of random sin¬ 
glet configurations that are completely unaffected by the 
cut: G = 1 and Sa = 0. The number of such config¬ 
urations _scales with a non-trivial exponent that can be 
deduced- 17 * 44 * from (|9| for example: P(n = 0, n e = 0) ^ 
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L ( 3 ^)/ 4 . This critical exponent can be rewritten as 
1 — ^ where cp = 1+ ^ is the golden ratio. The golden 
ratio appears in many different quantities in random sin¬ 
glet critical points, all of which being related in one way 
or another to the probability that a given bond be not 
decimated in the course of the RG. This yields 

M a ) 2* — ~r~' ( 18 ) 

a—^oo 4 

For a 1 , we therefore find G a ~ which is 

dramatically different from the very small typical scaling 
Gf yp ~ L~ a ^/ 4 . In this limit, the average G a is thus 
dominated by rare configurations with G = 1 that are 
unaffected by the cut. 

Of course, in general G a ~ A(a)L ^°^ where A(a) 
is a non-universal constant that will vanish in the limit 
ol — y oo so that strictly speaking, liuio^oo G a = 0 . How¬ 
ever, as we will see below, G a is dominated by rare con¬ 
figurations (and thus decays with an exponent ~ 3 ~^ ) 
already for a > 10. 


B. Multifractal Anderson Orthogonality exponents 
and generating functions 


We now turn to the calculation of the full multi¬ 
fractal spectrum n(a). In order to do so, we need 
to compute the full generating function P(u,v) = 
ne -P( w ^eKff e of the joint distribution P(n,n e ), 

with G a = T (t^, 2-^/2 ) • This is obviously a very hard 
problem, and we expect the calculation below to be only 
approximate - whereas our results up to this point should 
be exact. 

The generating function G(v) = F(l,v) = P ne v Ue 
is related to the entanglement spectrum and was com¬ 
puted (in an approximate way) in Ref. T 7 | The key idea 
was to assume that this generating function satisfies the 
following renewal equation as a function of RG time /jl 

POO P[1 

Gn(v)= dn' f(n') + v ( 19 ) 

Jfi JO 

where the singlet formation rate /(/i) is given by ©• 
This equation can be readily solved by Laplace transform, 
and one finds 


G(v) 


( 20 ) 


where we have used the fact that In T ~ i In L. In par¬ 
ticular, one can check that n e = Q'(y = 1 ) = | In L. 
We emphasize that the renewal equation ( 19 ) is ignoring 
memory effects beyond first order, multiple decimations 
etc. It is actually possible to show numerically that ( 20 ) 


does not hold exactly (see below), even if it provides a 
very good approximation of the exact result. 

The calculation of the generating function F(u,v) 
is even more intricate. To do this, we assume that 


the joint probability distribution P(n, n e ) factorizes as 
P(n, n e ) « P(n e )H(n — n e ), where H{n — n e ) is the prob¬ 
ability distribution for the number of singlets modified 
by the cut that were not crossing the cut in |T^). Those 
events are governed by the rate function g(/i) = e - ^ in 
eq. and should be described by a Poisson process 

since they are essentially independent. For a configu¬ 
ration with n e singlets crossing the cut, there are n e /2 
“active regions” where n — n e can increase, and we have 
seen above that typically, each one of these regions has 
n = 4 singlets modified by the cut, so that n — n e = 2 
for each active regions. We thus take H{n — n e ) to be a 
Poisson process with parameter n e /2 x 2 = n e , so that 
P(n, n e ) ~ P(n e ) ^_ n e y n ^~ ne . We can now evaluate the 
generating function F(u,v) as 

P(u,v) « ^P(n e ) (m;e w-1 ) ne = £(m;e w-1 ), (21) 

n e 

so that the full generating function decreases as a power 
law as 


P(u , v) 


L 


3 —\/5 + 4 u ve( u 

4 



( 22 ) 


One can check that this formula gives d u P(u , v)\ u V=1 = 
n = | In L and d v F{u,v )\ uv=1 = n e = | In L as it 
should. Going back to the multifractal Anderson orthog¬ 
onality spectrum, this gives 


3 - \/ 5 + 4 (Tr ex P( 2 a ~ 1 ) 

^ - «± - . (23 ) 


This formula is compatible with the exact limits (|T 
and (Il8|). 


V. NUMERICAL RESULTS 

In this section, we check our results numerically con¬ 
firming that random singlet critical points suffer from an 
Anderson orthogonality catastrophe. Admitting that the 
RSRG provides an accurate description of the ground- 
state of the random chains © and Q, the universal scal¬ 
ing of the overlap G can be very efficiently computed by 
implementing the RSRG procedure numerically. Starting 
from a given disorder configuration, we identify and deci¬ 
mate the strongest bonds by forming singlets. We repeat 
this procedure until we spanned the whole chain, and we 
end up with the groundstate wave function consisting of 
a collection of singlets. We compute the groundstates 
with and without the cut for each disorder realization, 
and we compute the overlap G as explained in Sec. |II| 
We draw the random couplings Ji E [ 0 , 1 ] from the dis¬ 
tribution P(J) = pp ji-i/w with W = 5 . Results for 
different system sizes L are shown in Fig. [ 3 a| averaged 
over more than 2 x 10 6 disorder realizations. We find that 
G ty P , G and G 20 clearly vanish algebraically with L, and 
the corresponding critical exponents are found to be in 



















FIG. 3: (a) Power law dependence of Gt yp and G a as a function of system size L. The data points obtained from a numerical 
implementation of the renormalization group agree very well with our analytic predictions (solid lines). Inset: Numerical check 
of this power-law dependence for a random-bond XX chain using free fermion techniques (solid lines: RSRG predictions, dashed 
line: overlap for a clean chain), (b) Multifractal Orthogonality spectrum G a ~ L ~^°The numerically extracted exponent 
n(a) agrees quite well with eq. (23) (blue line). Inset: G a against L for various values of a = 10 -4 ,..., 10 2 . 


very good agreement with ( 12 ) and ( 23 ). In particular, 
G 20 decreases with L very slowly with an exponent close 
to ( 18 ). By fitting G a as a function of L, we also mea¬ 


sured the exponent /i(a) for various values of a (Fig. [ 3 b|) . 
The numerical results are in relatively good agreement 
with the (approximate) formula (J 23 |, but although the 
finite size corrections seem to tend in the right direction, 
eq. ( [ 23 ] ) is clearly slightly off for a ~ 1 . We find similarly 
that eq. ( |20| ) for the entanglement spectrurrp] is ruled 
out by numerical results, despite being very close to the 
exact solution. We leave the possibility of deriving exact 
formula for these quantities for future work. 


Our results can also be checked numerically for the 
random XX chain (A = 0 in eq. Q) that can be mapped 
onto free fermions after a Jordan-Wigner transformation. 
Since the Jordan-Wigner transformation is nonlocal, in 
principle one has to be careful when comparing real space 
properties. However, the key quantities we study, such as 
entanglement and overlaps, are independent of whether 
one chooses the spin or fermionic representation; at most 
there are constant factors related to the boundary condi¬ 
tion used in defining the Jordan-Wigner string, but these 
do not modify the scaling behavior that is the focus of 
our study. We point out that our results apply equally 
well to the fermionic problem. For a free fermion Hamil¬ 
tonian H = J 2 ij hij c \ c ji the groundstate wavefunction 
can be obtained by diagonalizing the matrix hij and by 
filling the Fermi sea. The overlap G between two differ¬ 
ent groundstates can then be expressed as a determinant 
using Wick’s theorem. Unfortunately, we find that the 
numerical evaluation of such determinants becomes un¬ 
stable at strong disorder and/or large distances - that 
is, in the random singlet regime we are interested in. 


We note that even though similar free fermions meth¬ 
ods were used to check the scaling of the entanglement 
entropy in the random singlet regime^, the computa¬ 
tion of the entanglement entropy is also plagued by nu¬ 
merical instabilities at strong disorder^. Groundstate 
overlaps seem unfortunately even more sensitive to such 
instabilities, making it hard to access the random singlet 
physics. Nevertheless, by restricting our numerics to a 
regime (corresponding to moderate system size L) where 
the instabilities are negligible, we find a good agreement 
between numerical results for couplings drawn from a 
uniform distribution Ji G [0,1] and our analytic predic¬ 
tions (see inset in Fig. 3 a). 


VI. GENERIC PERTURBATIONS AND LOCAL 
QUANTUM QUENCHES 

In this section, we go back to more general arbitrary 
local perturbations and conjecture the form of the dy¬ 
namical response of a random singlet critical point to a 
sudden local perturbation (local quantum quench). 

We mostly focused above on a specific type of local 
perturbation that was amenable to analytic calculations: 
namely, a cut in the system. It is also instructive to 
consider more generic local perturbations. For example, 
we also computed (by implementing the RSRG proce¬ 
dure numerically) the overlap G = between 

the groundstate ^ of a disordered spin chain of size 
L = 2 N with random couplings Ji,..., Jjv? • • • Jl- i, and 
the groundstate of the same disorder configuration 
with a weakened central bond = XJn with A < 1 (the 
“cut” case considered in most of this paper then corre- 
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L/(l nA ” 1 ) 2 

FIG. 4: Crossover physics of the overlap G(L) for a local 
perturbation corresponding to a weakened link Jn = A Jn 
with A < 1. Once plotted against L/L* ~ L/ (in A x ) 2 , we 
find that the data for the overlap collapse up to a rescaling 
by a non-universal prefactor /(A) that we adjust freely. Inset: 
Universal collapse of the logarithmic derivative (no free 

parameter). 


global non-universal prefactor /(A) that we adjust freely) 
for different values of A (see Fig. [ 4 ^. To get rid of this 
non-universal A-dependent part in G(L), we also checked 
that the logarithmic derivative collapses onto a sin¬ 
gle universal curve without any free parameter onceplot- 
ted against L/L* ~ L/ (in A -1 ) 2 (see inset in Fig. |i). 


B. Local quantum quenches 


Let us also briefly comment on the generalization of our 
results to out-of-equilibrium setups. Let be the 

(sorted) eigenstates of the Hamiltonian Ha of a random 
spin chain and let ^ be the eigenstates of the Hamil¬ 
tonian Hb = Ha + V of the same system with a local 




(o) 


, we 


perturbation V. Starting from the initial state 

imagine a quench protocol where the local perturbation is 
suddenly turned on at time t = 0, so the wave-function 

of the system at time t reads |T(£)) = e 


— p—iHst 




(o) 


The time-dependent analog of the groundstate overlap 


g = <U 0) 


■v) 


sponding to A = 0 ). We find that the system suffers from 
an Anderson orthogonality catastrophe only in the “cut” 
case A = 0 . In other words, the nature of this algebraic 
decay of the groundstate overlap seems to be very differ¬ 
ent from the clean case where any arbitrary weak local 
perturbation is enough to induce an Anderson orthogo¬ 
nality catastrophe. For random singlet systems, we find 
that the overlap G remains finite in the thermodynamic 
limit L —»• oc for any value of A ^ 0 . 


A. Crossover for generic local perturbations 

Using the scaling of the random singlet fixed point, we 
introduce a characteristic length scale L\ (in A 1 ) 1 ^’ 
with ^> —\ associated with the energy scale of the weak¬ 
ened link. At high energies Q Jn (corresponding to 
system sizes L <C L\ ), the weakened bond can be es¬ 
sentially considered as Jn ~ 0, and our results for the 
cut should apply in this regime with the overlap G a scal¬ 
ing as L~^ a \ This is very intuitive physically: if L is 
small enough (or equivalently if A is small enough) so 
that L < La, the central bond with strength Jn will al¬ 
most surely be the last bond decimated during the RSRG 
procedure so there is no difference between a weakened 
bond and an actual cut. On the other hand, as the cut¬ 
off Ll is lowered and becomes of order Jn (L rsj L x ), the 
fact that the weakened bond is in fact non-zero becomes 
important, and the overlap G goes to a (non-universal) 
constant that depends on A as L L\. We verify 
this crossover scaling by collapsing G(L) curves (up to a 


G(t) = (U 0) \e iHAt e ~ iHBt | <U 0) ) , ( 24 ) 


which up to a phase is nothing but the overlap of the 
wavefunction at time t with the initial wavefunction at 
time t = 0 . We will consider the case where the lo¬ 
cal quench corresponds to either two random spin chains 
suddenly glued together, or to a single random spin chain 
suddenly cut in half at time t = 0. 

For pure gapless spin chains, the Loschmidt echo after 
such a local quench decayP^ as a power-law G(t) ~ 
where r] is the same exponent as the orthogonality expo¬ 
nent characterizing the decay of the “static” groundstate 
overlap G ~ L~ r] as a function of the system size 39 - this 
can essentially be traced back to the symmetry between 
time and space (dynamical exponent z = 1) of these sys¬ 
tems. The Loschmidt echo seems quite hard to compute 
using RSRG but it is natural to expect the behavior of 
G(t) in the time dependent setup to be related to the 
static scaling G a ~ us i n g the low-energy relation 

In t VL (see eq. © with t 1 = E acting as an energy 
scale). We therefore conjecture that for local quantum 
quenches, the long-time behavior of the Loschmidt echo 
should be given by 


G(t) a 


1 

(lnt) 2 ^’ 


( 25 ) 


where fi(a) is the same exponent as the one appearing 
in the static groundstate overlap G a ~ L~^ a \ Using 
the same argument for the entanglement entropy yields 
a very slow growth S(t) ~ lnlnt after a local quench, 
consistent with Ref. [Ml Note that the behavior after a 
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FIG. 5: Slow decay of the Loschmidt echo after two random 
XX chains of size L/2 are suddenly connected by a random 
bond (local quench). 


to zero by convention - this can always be done by shift¬ 
ing the frequencies uj — ujq u. For our local quench 
consisting of two random singlet spins chains glued to¬ 
gether at time t = 0 , the average Loschmidt echo is given 
by eq. ( [ 25 ] ) with a = 1 so that the work distribution 
(absorption spectrum) scales at low frequencies as 


A(uj) ^ 




( 27 ) 


with r] = /i( 1 ). Note that this singularity is integrable, 
compatible with the fact that f dcjA(cu) = 1 is finite. 
The very slow decay ( 25 ) of the Loschmidt echo thus 
translates into a highly singular (barely integrable) be¬ 
havior at low frequencies. It would be interesting if the 
methods of Refs. 1151521 to compute equilibrium dynam¬ 
ical response functions in random spin chains could be 
extended to recover this result. 


global quench would be dramatically different: consider¬ 
ing the entanglement entropy for instance, the entangle¬ 
ment growth starting from a product state for the ran¬ 
dom Ising chain becomes S(t) ~ (lnt) Q in the presence 
of interactions breaking integrability 22 . 

We checked the scaling ( [ 25 ] ) numerically in the ran¬ 
dom XX chain (with couplings drawn from the uniform 
distribution on [0,1]) where the local quench consists in 
suddenly connecting two chains of size TV by a random 
coupling drawn from the same distribution (Fig. §• By 
computing the Loschmidt echo as a determinant using 
Wick’s theorem, we found a long time behavior of G(t) a 
compatible with eq. ( 25 ), with the exponent fi(a) consis¬ 
tent with our results in the static case. 


C. Work distribution 


Another quantity of interest is the Fourier transform 
A{oo) of the Loschmidt echo ( [ 24 ] ) G(t) = f duoe~ luJt A(u) 
that corresponds physically to the distribution of the 
work done during the quantum quench 47 51 




A(u) = ^ 




‘s (u - E[ n) + ^ 0) ) . ( 26 ) 


In optical absorption experiments realizing quantum 
quenches, A(w) corresponds to the absorption spectrum 
where the energy hu of a photon naturally coincides with 
the work needed to perform the quantum quenchPSEHl 
In clean systems, the power-law decay of the Loschmidt 
echo with time G(t) ~ then translates into an edge 
singularity A(u) ~ ^(cj )^ -1 similar to the X-ray singu¬ 
larity in metalsT 2 ', with 0{pc) the Heaviside step function. 
Notice that in this expression, we set the minimal work 
required to perform the quantum quench ujq = e[°^ —E^ 


VII. CONCLUSION 


We have studied analytically and numerically the re¬ 
sponse of random singlet quantum critical points to lo¬ 
cal perturbations. We showed that such systems suffer 
from an Anderson orthogonality catastrophe very sim¬ 
ilar to that in metals, with the overlap of the ground- 
states with and without a cut in the system decaying 
algebraically with system size. We showed that the or¬ 
thogonality catastrophe is multifractal in the sense that 
different powers of the disorder-averaged overlap G scale 
with non-trivial powers of system size. We also discussed 
the implications for more general local perturbations and 
for the nonequilibrium dynamics of such systems after a 
local quantum quench. 

Using real space renormalization group techniques, we 
were able to compute the corresponding orthogonality 
exponents exactly in some limits, and approximately for 
the full multifractal spectrum. Whether exact expres¬ 
sions can be obtained from the RSRG for the multifractal 
orthogonality spectrum and for the entanglement spec¬ 
trum 17 remains an open question. We also note that the 
nature of the orthogonality catastrophe discussed here 
is quite different from that in metals: in particular, we 
showed that the overlap remains finite (on averag e) for 
more generic, “softer” local perturbations (see Sec. VIA 
above). Similarly, coupling a random-singlet system to 
a quantum impurity will not lead to an orthogonality 
catastrophe, contrary to the metallic case. The physical 
mechanisms underlying the orthogonality catastrophe in 
random-singlet systems are also quite different from those 
at play in metals: within the RSRG, a physical cut in 
a random-singlet system leads to an avalanche of local 
rearrangements of singlets that in turn leads to the or¬ 
thogonality catastrophe. It would be interesting to study 
this orthogonality catastrophe in the free fermions case 
(XX chain) using different approached^, and to perhaps 
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relate it to the more usual language of phase shifts or to 
the singular form of the density of states at energy 5 = 0. 
It would also be very interesting to investigate whether 
the orthogonality catastrophe discussed here generalizes 
to highly excited states - in the context of many-body 
localization - for systems that can be tackled using sim¬ 
ilar RSRG techniqueJ^, and if it can be related to the 
discussion in Ref. EU 
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Note that this can be checked independently. The rate § 
leading to the RG duration In jA = 3 relies on two ingre¬ 
dients: (a) a flow equation for the central coupling similar 
to (10) with a factor 2 in front of the second term to ac¬ 
count for the decimations of its two neighboring bonds, (b) 
the initial condition at scale To at which the central bond 
was just decimated and is therefore strongly suppressed. At 
the scale T a > To however, it is natural to consider that 
the central bond is generic and is therefore distributed ac¬ 
cording to Qr a (fl) = p-e _/3//r<x . Solving the flow equation 
starting from this initial condition at T = T a , we find that 


In = 2 consistent with In ^ = 1. 
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